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Abstract — Under the deregulation of generation market in China, 
all distributed generators will particular in electric power 
bidding. Therefore power purchase cost optimization (PPCO) 
problem has been getting more attention of power grid Company. 
However, under the competition principle, they can purchase 
power from several of power plants, therefor, there exist 
continuous and integral variables in purchase cost model, which 
is difficult to solve by classical linear optimization method. An 
improved differential evolution algorithm is proposed and 
employed to solve the PPCO problem, which targets on minimum 
purchase cost, considering the supply and demand balance, 
generation and transfer capability as constraints. It yields the 
global optimum solution of the PPCO problem. The numerical 
results show that the proposed algorithm can solve the PPCO 
problem and saves the costs of power purchase. It has a widely 
practical value of application. 
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I. Introduction 

With the implementation of separating the power plants 
from grid company and generation bidding in power market in 
China, it is important for the grid company to optimize the 
purchasing schedule to reduce the cost of grid company[l-3]. 

According to the different settlement rules, the objective 
function of generation-side electricity market is divided to two 
classes. One is to minimum the cost based on the unified 
marginal settlement, while the other is minimum cost based on 
the practical generators bidding[4,5]. The second power 
purchase cost optimization(PPCO) problem is studied in this 
paper. The PPCO problem is studied for some years. DUAN 
established a purchase optimization model considering the 
generation capability and line loss, and solved this model by 
penalty function [1]. But he only discussed the PPCO problem 
on the protection and coordination principle. CHEN et al. 
established purchase optimization model on the market 
principle and solved it by Lagrange relaxation method [4]. 
TAN et al. discussed the multilevel purchase optimization 
model of long-term, mid-term and short-term electricity market 
and analyzed the relationship of them, but they did not give the 
solution method [6]. ZHANG et al. established a purchase 
optimization model with the constraint of supply and demand 


balance and constraint of generation capability. They solved 
this model by the PSO but they did not consider the transfer 
capability constraint and did not satisfy the constraint of supply 
and demand balance [2]. ZHANG et al. established a purchase 
optimization model with the constraints of supply and demand 
balance, generation capability and transfer capability. They 
solved it with PSO and obtain a satisfied solution [3]. 

Differential Evolution (DE) algorithm [7, 8] is a heuristic 
global search method which has been employed in many fields 
such as chemical industry, machine design, information 
processing, biological information, geophysical 
inversion[9,10,ll]. The DE algorithm has employed to solve 
the PPCO problem recently, but did not consider the 
constraints of transfer capability and supply and demand 
balance [12]. 

The purchase optimization model is established with the 
constraint of the supply and demand balance, generation 
capability, and transfer capability firstly. And an improved DE 
algorithm with the dynastic mutation factor and crossover 
factor is proposed secondly. Then the steps of algorithm for the 
PPCO problem based on the improved DE algorithm thirdly, 
and the correctness of the proposed algorithm is verified by 
numerical experiment. 

II. Power Purchase Model 

The goal of power purchase optimization is to minimum the 
cost of purchasing particular quantity of electricity considering 
some constraints. 

The power purchase cost optimization problem can be 
modeled as follow: 

n 

MinF, F = Y J C i P i 

i = 1 

S.T. £d -A. m=P d 

i = 1 

P ' = 0 or P < P < P 

i i min — i — i max 

- L < L < L 

i max i i max 

( 1 ) 
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where C i ,P i ,AP i is the electricity price, the purchase 
quantity of electricity and the line loss of i-th power plant 
respectively; P d is the particular quantity of electricity should 

be purchased; P inAn and P irmx is the minimum and maximum 

generation capacity of z-th power plant respectively; L /max is 

the transmission capacity of the line from z-th power plant to 
grid company; F is the objective function of power purchase 
cost optimization. 

III. Improved DE algorithm 

A. Basic DE algorithm 

DE algorithm first initialize is a parallel direct search 
method which utilizes NP n-dimensional parameter vectors 

,x? ,...,x° p J as a population for each generation. The 

initial vector population is chosen randomly and should cover 
the entire parameter space. Then DE algorithm update 
population by the operation of mutation, crossover and 
selection. 

(1) Mutation 

DE generates new parameter vectors by the operation 
named mutation which add the weighted difference between 
two population vectors to a third vector. 

vr=x k ri + F(x k r2 -x k r3 ) ( 2 ) 

where x k = (x^,xf 2 ,...,x^)is the z-th vector of &-th generation, 
r 19 r 2 ,r 3 e {1,2,..., NP} , r 1 Ar 2 ^r 3 ^z, F > 0 is the mutation 
factor. 

(2) Crossover 

In order to increase the diversity of the perturbed parameter 
vectors, crossover is introduced as follow. 

k+l [ vj +1 , rand ( ij ) < CR or j = rand (j); 

u it ={ . (3) 

[ Xy , rand (ij) > CR or j ^ rand ( j ); 


If the algorithm convergence or the max iteration is 
achieved, algorithm finish, otherwise, begin the next iteration 
of mutation, crossover and selection. 

B. Improved DE Algorithm 

The mutation factor E and crossover factor CR is critical to 
the performance of basic DE algorithm. The large F and CR is 
beneficial to the global search ability, but it is harmful to the 
local search ability leading to the poor convergency speed, and 
vice versa. So we improve the DE algorithm in the following 
ways. 

(1) Adaptive mutation factor 

The following adaptive mutation factor is introduced in 
order to improve the global search ability at the early stage and 
the local search ability at the later stage of the algorithm. 

t (P -p \ 
rg r v max min z 


(5) 

(2) Increasing crossover factor 

DENG and RAN [13] propose the following strategy of 
increasing crossover factor to achieve the balance between 
global search ability at the early stage and local search ability 
in the later stage of the algorithm. 

CR = (CR. - CR um )[(— ) 2 - 2 — + 1] + 0? 

' min max ^ l \ ^ max 

max ± max 

( 6 ) 

(3) Random mutation 

In order to avoid the stagnation of DE algorithm at its last 
stage, a random mutation method is introduce. If an individual 
are stagnated in certain iterations, a random mutation is happen 
in this individual, that is random initializing a individual to 
replace the stagnated one. The random mutation method is 
described as follow: 

if F(xj) = F(x\ +l ) = F(x\ +2 ) = ■■■ = F(x\ +P )and F(x\) * F* 


where rand(ij ) is a uniform random number generator 
with outcome rand(ij) e [0,1] ; CR is crossover factor 
CR g [0,1] which has to be determined by the user; rand(j) 
is a randomly chosen index rand(j) e 1,2,...,£) , which 

k~\~\ k~\~\ 

ensures that U- gets at least one parameter from . 

(3) Selection 

To decide whether or not the U^ +l should become a 

k 

member of generation k + 1 , it is compared to the X tj using 
the greedy criterion as follow. 

r«* +1 ,/(M* +1 ) </(**); 

‘ ' ‘ (4) 

' I x , Otherwise. 


then x\ +p+x = x^ B + rand (0, 1) * (x max - ) 

(7) 

Where F is fitness of global best individual, p is max 
stagnation iteration allowed, (x min , A max ) is the random search 
space. 

We define the improved DE algorithm, which employ the 
adaptive mutation factor and increasing crossover factor 
strategy, named as DEI. And we define the improved DE 
algorithm, which employ the adaptive mutation factor, 
increasing crossover factor, and random mutation strategy, 
named as DE2. 
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IV. Solving the PPCO problem with Improved DE 

ALGORITHM 


A. Strategy for the constraints 


The generation and transfer capacity constraints are 
transformed to boundaries of parameter vector. The balance of 
supply and demand constraints is considered using the penalty 
function method and the PPCO problem is transformed from a 
optimization problem with constraints to a non-constraint 
optimization problem. The objective function is following: 


MinF , F^C.P.+a 

i = 1 


£fl- A p,yp-p^„ 

_ i = 1 


( 8 ) 


where P Qxpect is the particular power should be purchased; 


a is penalty factor; the variables be optimized P t , is the 
power purchased from z-th power plant. 


B. Algorithm steps of improved DEfor PPCO 

The algorithm steps of improved DE for PPCO problem are 
following: 

Step 1: given the upper and lower boundaries of the variables 
be optimized P t , the population NP, the upper and lower 
boundaries of mutation factor F and crossover factor CR; 


Step 2: initialize the population of parameter vector randomly 
according to the upper and lower boundaries of P t ; 


Step 3: calculate the objective function value of P t , according 
to formula (8); 

Step 4: execute the mutation, crossover and selection according 
to formula (2), (3), (4) respectively, and obtain the next 

generation population [x^ +1 , x \ +1 , . . . , X k ^p ] ; 

Step 5: update the mutation factor and crossover factor 
according to the formula (5), (6); 

Step 6: if the algorithm stagnate certain iterations, random 
mutation a individual according to the formula (7); 


Step 7: calculate the objective function value of P { , according 
to formula (8); 

Step 8: if the convergency and the max iteration is not 
achieved, go to step 4, begin next iteration; 

Step 9: output the solution, finish. 

V. Numerical Result and Analysis 


A. Benchmark Results 

We use following five benchmark functions to test our 
improved DE algorithm. 

(1) Sphere 

/ i (x) = Jxf, M = 30,1^1 < 100 (9) 

1=1 

(2) Generalized Schwefel’s 

/ 2 w=E (a. sin(^jx^j)) , n = 30,|x| < 500 (10) 

i = 1 

(3) Rastrigr 


/ 3 (x) = -lOcos^;^.) +10] ,77 = 30, |x.| <32 (11) 


(4)Griewank 

/ 4 M = ~ H cos( 4) = 1 ’ ” = 30,|x I ^ 600 (12) 

,=1 V/' 


4000 tt 


(5)Ackly 


(i 

\n Ti 


1 ^ 


f 5 (x ) = -20exp[-0.2 — ] _ex P — ^cos(2^x.) 


i=i 


+ 20 + e 


y 


,n = 30,|jc.| <5.12 (13) 

The parameters of standard DE algorithm are set as follows: 
CR = 0 . 9 , F = 0.5 . The Parameters of improved DE 

algorithm are set as follows: C R mn = 0 . 1 , CR mdx = 0 . 9 ? 


F m - n = 0 . 3 , F nrAx = 1.2 . Both algorithms’ number of 

individuals are 60, max iteration are 2000. The test results are 
listed in the table 1. 


TABLE 1. Benchmark on DE and improved DE algorithm 




f i 

u 

u 

h 

fs 

DE 

Average 

6.47e-31 

-9.74e+3 

102.02 

7.40e-4 

5.51e-15 

Stand, dev. 

1.33e-39 

900.94 

35.19 

0.0022 

1.63e-15 

Best 

1.51e-32 

-1.12e+4 

10.94 

0 

4.45e-15 

Worst 

7.34e-30 

-7.43e+3 

183.66 

0.0074 

7.99e-15 

Runtime 

85.92 

131.89 

107.17 

118.06 

97.97 

Improved 

DE(DEl) 

Average 

4.98e-5 

-1.26e+4 

3.32e-5 

1.77e04 

0.0015 

Stand, dev. 

9.42e-5 

2.30e+4 

5.59e-5 

2.29e-4 

0.0014 

Best 

7.66e-8 

-1.26e+4 

8.63e-9 

8.39e-8 

8.85e-5 

Worst 

2.43e-4 

-1.26e+4 

3.08e-4 

8.67e-4 

0.0055 

Runtime 

262.24 

260.94 

278.27 

279.52 

284.47 

Improved 

DE(DE2) 

Average 

4.13e-5 

-1.26e+4 

1.67e-5 

8.23e-05 

0.0012 

Stand, dev. 

6.81e-5 

12.4e-4 

6.95e-5 

1.34e-4 

0.0013 
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Best 

9.61e-12 

-1.26e+4 

7.59e-10 

2.02e-8 

2.14e-5 

Worst 

4.72e-4 

-1.26e+4 

3.01e-4 

4.89e-4 

0.0050 

Runtime 

290.94 

274.67 

279.70 

323.17 

267.74 


The results listed in table 1 indicate that the both DEI and 
DE2 can improve the accuracy of fitness of benchmark 
functions, while DE2 is better than DEI at small degree. 

B. PPCO Results 

The problem of grid company purchases 200GWh power 
from 5 power plant is considered. The model of power supply 
is illustrated as Figure 1. 

There is a backup line between plant 2 and 3, which is 
opened normally. When the line between plant 2 and grid 
company or between plant 3 to grid is broken, the backup line 
is closed and the plant 2 and plant 3 share one line to transfer 
the power. The line loss of the backup line between plant 2 
and plant 3 is 0.002. 

The emulation program is written using Matlab m 
language. The parameters are chosen as following: 
D=5, NP =40, F_ = 0 . 3 , = 1.2 , CR mn = 0 . 1 , CR^ = 0.9 , max 
iteration is 5000. The price, line loss, generation and transfer 
capability are listed in table 2. 

Under the normal condition, that means the transfer 
capability is larger than the generation capability. The results 
of purchase optimization under the protection and 
coordination principle are listed in the table 3, compared with 
the results of literature [2] [3] [4]. The results show that the 
cost of our algorithm is same as that of literature [3], but 
larger than those of the literature [2] [3] slightly, because the 
results of ours and literature [3] is obtained under the 
constraint of supply and demand balance restrictively whereas 
those of literature [2] [4] is obtained with some deviation of it. 

The results of purchase optimization under the normal 
condition and marketing principle are listed in the table 4, 
compared with the results of literature [2] [3] [4]. The results 
show that the cost of our algorithm is same as those of 


literature [3] [4] which are satisfied with the constraint of 
supply and demand balance restrictively, while the result of 
literature [2] is not satisfied that constraint restrictively. 

When the line between plant 3 and grid company is broken, 
the backup line between plant 2 and 3 is close and the two 
plants share one transfer line between plant 2 and grid 
company. The results of purchase optimization under such the 
abnormal condition are listed in the table 5. The optimization 
results are same as those of literature [3]. Under the protection 
and coordination principle is 27.6789 million Yuan(Y), and 
under the marketing principle is 27.3541 million Yuan(Y) . 

The numerical results show that the proposed improved 
DE algorithm has the advantage of fine global optimize 
ability. It can be employed to solve the PPCO problem and 
reduce the purchase cost obviously. 



Fig. 1. power supply model of 5 power plants 


TABLE 2. Parameters of powe r supply model of 5 power plants 



c,(¥/kWh) 

AP, 

Pi^JGWh 

P imax /GWh 

L^JGWh 

plant 1 

0.10 

0.0882 

43.2 

86.4 

100 

plant 2 

0.12 

0.0722 

21.6 

64.8 

90 

plant 3 

0.15 

0.0451 

21.6 

43.2 

60 

plant 4 

0.18 

0.0422 

14.4 

43.2 

60 

plant 5 

0.20 

0.0554 

14.4 

28.8 

40 


TABLE 3. The results of power purchase optimization under the protection and coordination principle 



literature [2] 

literature[3] 

literature [4] 

our results 

plant 1 

86.4000 

86.4000 

86.3979 

86.4000 

plant 2 

64.8000 

64.8000 

64.6384 

64.8000 

plant 3 

35.2963 

35.6356 

35.2467 

35.6356 

plant 4 

14.4000 

14.4000 

14.5145 

14.4000 

plant 5 

14.4000 

14.4000 

14.4549 

14.4000 

cost/million Y 

27.1825 

27.2333 

27.19498 

27.2333 

supply/demand balance 

199.6760 

200.0000 

199.6391 

200.0000 
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TABLE 4. The results of power purchase optimization under the marketing principle 



literature[2] 

literature [3] 

literature[4] 

our results 

plant 1 

86.4000 

86.4000 

86.4000 

86.4000 

plant 2 

64.8000 

64.8000 

64.8000 

64.8000 

plant 3 

43.2000 

43.2000 

43.2000 

43.2000 

plant 4 

20.7218 

21.0601 

21.0601 

21.0601 

plant 5 

0.0000 

0.0000 

0.0000 

0.0000 

cost/million Y 

26.6269 

26.6868 

26.6868 

26.6868 

Supply/demand 

balance 

199.6760 

200.0000 

200.0000 

200.0000 


TABLE 5. The results of power purchase optimization under the abnormal condition 



literature [3] 

our results 

protection 
and coordination 
principle 

marketing 

principle 

protection and 
coordination 
principle 

marketing 

principle 

plant 1 

86.4000 

86.4000 

86.4000 

86.4000 

plant 2 

64.8000 

64.8000 

64.8000 

64.8000 

plant 3 

25.2000 

25.2000 

25.2000 

25.2000 

plant 4 

25.5659 

39.7675 

25.5659 

39.7675 

plant 5 

14.4000 

0.0000 

14.4000 

0.0000 

cost/million Y 

27.6779 

27.3541 

27.6779 

27.3541 

Supply/demand balance 

200.0000 

200.0000 

200.0000 

200.0000 


VI. Conclusion 

An improved DE algorithm, which introduce new update 
method of mutation and crossover factor to improve the 
balance between the fine global search ability at the early 
stage and the fine local search ability at the later stage, is 
employed to solve the power purchase cost optimization 
(PPCO) problem with generation and transfer capability 
constraints. The results of five plants power supply model 
show that the proposed algorithm can find the global optimum 
of this PPCO problem and reduce the purchase cost obviously. 
It has a widely practical value of application. 
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